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Abstract. This paper develops three polynomial-time pricing techniques for Eu- 
ropean Asian options with provably small errors, where the stock prices follow 
binomial trees or trees of higher-degree. The first technique is the first known Monte 
Carlo algorithm with analytical error bounds suitable for pricing single-stock options 
with meaningful confidence and speed. The second technique is a general recursive 
bucketing-based scheme that can use the Aingworth-Motwani-Oldham aggregation 
algorithm, Monte-Carlo simulation and possibly others as the base-case subroutine. 
This scheme enables robust trade-offs between accuracy and time over subtrees of 
different sizes. For long-term options or high-frequency price averaging, it can price 
single-stock options with smaller errors in less time than the base-case algorithms 
themselves. The third technique combines Fast Fourier Transform with bucketing- 
based schemes for pricing basket options. This technique takes polynomial time in the 
number of days and the number of stocks, and does not add any errors to those already 
incurred in the companion bucketing scheme. This technique assumes that the price 
of each underlying stock moves independently. 
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1. Introduction 

A call (respectively, put) option is a contract assigning its holder the right, but not 
the obhgation, to buy (respectively, sell) a security at some future time for a specified 
strike price X [|T3|. If the holder exercises her right, the other party in the contract, the 
writer, is obligated to assume the opposite side of the transaction. In exchange for this 
right, the holder pays the writer an option price P. The security in this contract can 
be any financial asset; for the purpose of this paper, we restrict it to a single stock or a 
portfolio of stocks. An option in the latter case is commonly called a basket option; for 
clarity, we call an option in the former case a single-stock option. 

Options are popular financial instruments for a variety of trading strategies. For 
example, options can be used to hedge risk. As protection from a potential price fall in 
a stock price, one can purchase a put on the stock, thereby locking in a minimum sell 
price. On the other hand, options can provide additional income for stockholders who 
write calls on their holdings; of course, this strategy carries the risk of being forced to 
sell the stock should the calls be exercised. 

An option is valid until its expiry date. For a European option, the holder may 
exercise it only on the expiry date. For an American option, the holder may exercise it 
on any date up to and including the expiry date. The payoff of an option is the amount 
of money its holder makes on the contract. A European call is worth exercising if and 
only if S" > X, where S is the stock price on the expiry date. The payoff of the call is 
{S — X)'^ = max(5' — X, 0). For an American call, S is set to the stock price at the 
exercise time. An Asian option comes in European and American flavors, depending on 
when it may be exercised. For an European Asian call, if A is the average stock price 
over the entire life of the contract up to the expiry date, the payoff is {A — X)^. For 
an American Asian call, A is set to the average stock price up to the exercise date. The 
payoffs of puts can be symmetrically defined. 

The price P of an option is the discounted expected value of the payoff 
with an appropriate martingale measure. Because of the popularity of options. 



pricing techniques for computing P have been extensively researched [y, |, |, |^,[T9 
Generally, it is more difficult to price a basket option than a single-stock option. To 
compute P, the price movement of each individual stock needs to be modeled. Typically, 
it is modeled as Brownian motion with drift. Using a stochastic differential equation, P 
can then be computed via a closed-form solution to the equation. When a closed-form 
solution is not known, various approaches are used to find an approximate solution. 
One class of approaches involves approximating the solution using numerical methods. 
Other approaches approximate the Brownian motion model with a discrete model, and 
use this model to approximate P. One such discrete model is the binomial tree model. 



due to Cox, Ross, and Rubenstein ; see Section |^ for the definition of the model. 

This paper develops three polynomial-time pricing techniques with provably small 
errors. The remaining discussion makes the following assumptions: 

(i) The option in question is an European Asian single-stock or basket call. 
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(ii) Our task is to price the call at the start of its contract life. 

(iii) The price of each underlying stock follows the binomial tree model. 

(iv) In the case of a basket call, the price of each underlying stock moves independently. 

Our results generalize easily for puts, for a later time point than the start of the contract 
life, and for trees with higher degrees than two. The cases of American options and of 
interdependent stocks remain open. 

Monte Carlo simulation has been commonly used in the financial community. 
Despite this popularity, most reported results on error bounds are experimental or 
heuristic Our first technique is the first known Monte Carlo algorithm that 

has analytical error bounds suitable for pricing a European Asian call with meaningful 
confidence and speed. As shown by Theorem ^, the number of simulations required is 
polynomial in (PI) the logarithm of the inverse of the error probability and (P2) the 
inverse of the price error relative to the strike price but is exponential in (El) the square 
root of the number of underlying stocks and (E2) the volatility of these stocks over the 
call's life. In particular, the algorithm is reasonably fast and accurate for a single-stock 
European Asian call with reasonable volatility. 

Monte Carlo simulation is a randomized technique, and thus there is always a 
nonzero probability that the price obtained by polynomial-time Monte Carlo simulation 
is not accurate enough. The aggregation algorithm of Aingworth, Motwani, and Oldham 
(AMO) []I| is the first polynomial-time algorithm for pricing single-stock European 
Asian calls and other path-dependent options with guaranteed worst-case price errors. 
The AMO algorithm is based on a simple yet powerful idea called bucketing. Our 
second technique is a general recursive bucketing-based scheme that can use the AMO 
algorithm, Monte-Carlo simulation, and possibly others as the base-case subroutine. 
This scheme enables robust trade-offs between accuracy and time over subtrees of 
different sizes. For long-term options or high-frequency price averaging, it can price 
single-stock European Asian calls with smaller error bounds in less time than the 
base-case algorithms themselves. In particular, as implied by Theorem |14|, given the 
same runtime, this recursive scheme prices more accurately than the AMO algorithm; 
similarly, given the same accuracy, the scheme runs faster than the AMO algorithm. 

This recursive scheme works for calls written on a single stock. Our third technique 
combines Fast Fourier Transform (FFT) and bucketing-based schemes to price 
basket calls and is applicable to European Asian calls as well as others. As shown 
in Theorem |18|, this technique takes polynomial time in the number of days and the 
number of stocks, and does not add any errors to those already incurred in the companion 
bucketing schemes. 

The remainder of this paper is organized as follows. Section |^ reviews the binomial 
tree model and basic definitions. Section ^ describes the new Monte Carlo algorithm. 
Section ^ details the recursive scheme. Section |^ gives the FFT-based technique for 
pricing basket calls. Section concludes the paper with directions for further research. 
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2. The Binomial Tree Model 

A binomial tree T is a recombinant binary tree. If n is tlie deptli of T, T lias t + 1 
nodes at depth t, for < t < n. For < z < t, let T[t, i] (or simply [t, i] if T is obvious 
from context) be the i-th node at level t of T. For t > 0, T[t,0] and T[t,t] have one 
parent each, T[t — 1,0] and T[t — l,t — 1] respectively. For < i < t, T[t,i] has two 
parents, T[t — l,i — 1] and T[t — The number of nodes in T is ^ 

Given a stock in the binomial tree model, the stock price is assumed to follow a 
geometric random walk through T. Time is divided into n equal periods, with the 
root T[0, 0] corresponding to time t = 0, when the option is priced, and the leaves 
T[n, ■ ] corresponding to time t = n, the expiry date of the option. Let s(T[t,i]) (or 
simply s(t,i)) be the stock price at node T[t,i]. At each time step, the stock price 
s{t,i) rises to sit + l,i + 1) = u ■ s{t,i) — an uptick — with probability p or falls to 
s{t + 1,^) = d ■ s{t,i) — a downtick — with probability q = 1 — p. Letting r denote the 
risk- free interest rate, the parameters u and d satisfy < d < l+r < u and are typically 
taken to be m = ^ = e""/^, where a is the n-period volatility, or standard deviation, of 



the stock price |13|. Although the probability p of an uptick is not known in general. 



for the purposes of pricing options, we can use the risk-neutral probability model | |13|| , 
which states that p = -^^^i^, where r is the risk-free interest rate for one period. This 
makes the expected return on the stock over one period, pu -\- (1 — p)d, equal to the 
risk-free return, l + r. 

Let VL be the sample space of paths u = {ui, . . . ,ujn) down T, where each 
ujt G { — 1, 1}, with —1 corresponding to a downtick and 1 corresponding to an uptick. 
Given uj ^ Q and < t < n, let T[t, uj] be the unique node at level t that uj passes 
through. Similar to the notation introduced above, we let s(T[t,uj]) = s(t,uj) be the 
price at node T[t,uj]. 

We define the random variables Yi, ... ,Yn on Q hj Yt{uj) = ut, the t-th component 
of UJ. We define the probability measure 11 on to be the unique measure for which the 
random variables Fi, . . . , F„ are independent, identically distributed (i.i.d.) with P(Yi = 
I) = p and P(Yi = —1) = q. We start with an initial fixed stock price Sq = s(T[0, 0]). 
For 1 < t < n, the stock price St is a random variable defined by St = S'o'U^*=i^v 
From the structure of the binomial tree, we have Pr(S't = s{t,i)) = (*)p*(l — py~\ 
where < i < t. The running total of the stock prices is defined by Tt = Yll=o 
and the running average is At = Tt/{t + 1). For G we let Stiuj) = SQU^i=^^^^'^\ 
TtH = Elo •S'.(^), and Aiuj) = Tt{uj)/{t + 1). 

Recall that X is the strike price of a European Asian call. Using the above notation, 
the price of this call is 

E((^„ - X)+) = ^E((T„ - (n + 1)X)+) = ^E( max(T„ - [n + l)X, 0)). 

For large n, it is not known if this quantity can be computed exactly because the 
stock price can follow exponentially many paths down the tree. Below, we show how to 
estimate E((T„ — (n + l)X)+), from which the price of the option can be easily computed. 
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3. A New Monte Carlo Algorithm 

Monte Carlo simulation methods for asset pricing were introduced to finance by Boyle 
They are very popular in pricing complex instruments, particularly path-dependent 
European-style options. These methods involve randomly sampling paths u & Q 
according to the distribution 11 and computing the payoff {An{uj) — X)+ on each sample. 
Suppose samples u^, . . . , uj^ are taken from Vt. The price estimate of the call is 

N 



^ N 

1=1 

The accuracy of this estimate depends on the number of simulations and the variance 
of the payoff: the error bound typically guaranteed by Monte Carlo methods is 
0{t/^/N). Generally is not known, and itself must be estimated to determine the 
error bound of fi. 

A number of techniques are used to reduce the error of fi. For example, use 
the control variate technique, which ties the price of the option to the price of another 
instrument (the control variate) for which an analytically tractable solution is known. 
The simulation then estimates the difference between the option price and the control 
variate, which can be determined with greater accuracy than the option price itself. The 
antithetic variate method |^, 12] follows not only randomly sampled paths u down the 



binomial tree, but also the "mirror images" of each u. None of these techniques have 
known analytical error bounds for fi. 

Below we use concentration of measure results in conjunction with Monte Carlo 
simulation to estimate the price of a European Asian call and derive analytical error 
bounds for this estimate. The error bounds are in terms of the strike price X of the 
option and the maximum volatility o"max of the underlying stocks. 

3.1. Analytical Error Bounds for the Single-stock Case 

Let C = 6^(1^^^"). In this section, we show that if {n + 1)X/C is "small", E((T„ - {n + 
1)X)+) is close to E(T„ - {n + 1)X) = E(T„) - (n + 1)X (Theorem |(1)). The option 
is deep-in-the-money and will probably be exercised. Since a closed-form formula exists 
for E(T„), as described in Lemma |^, E(T„ — {n + 1)X) can be computed exactly, and 
our algorithm uses it as our estimate for E((T„ — (n + 1)X)+). On the other hand, if 
(n + 1)X/C is not small, the variance of (T„ — (n + 1)X)+ can be bounded from above 
(Theorem |^(2)) and our algorithm estimates its expectation with bounded error using 
Monte Carlo simulation. 

We first give some theoretical results, then show how these results can be used in 
our Monte Carlo algorithm, BoundedMC. We begin by finding bounds A and B such 
that T„ G [A, B] with high probability. Our main tool for this is an inequality arising 
from the theory of martingales, usually known as Azuma's Inequality [@, We use 



the inequality in a form taken from Frieze |TT|]. Suppose we have a random variable 



U = U{V), where V = (Vi, . . . , K;) and, for 1 < i < 77, Vi is chosen independently from 
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the probability space fli. For V,W e Q = fli x ■ ■ ■ x Q^, write V ~ W if V and W 
differ in exactly one component. 

Theorem 1 (Azuma's Inequality) Suppose that V implies \U{V) — U{W)\ < 
a, for some constant a > 0. Then, for any real b > 0, 

Pr(|[/-E(f/)| > 6) < 2e ^. 

Lemma 2 Let C = e'^*^''^^"). For any A > 0, we have 

Pr (r„ < Ce-''^ or T„ > Ce''^) < 26"^'/^ 
where a is the volatility of the stock. 

Proof: We apply Azuma's Inequality to the random variable U = lnT„. Let 
ujx,ujy & Q he two paths that differ in a single component, say on day t. Assume 
without loss of generality that ux has an uptick on day t while uy has a downtick. 
Then < and I In ^4^1 = \\nTJuJx) - InTJooyM < 2\nu = 2^. 

Applying Azuma's Inequality with a = 2^ and b = aX, we have 

Pr(| lnT„ - E(lnT„)| > aX) < 2e"^'/l 
The claimed result follows from the fact that 

I lnT„ - E(lnT„)| > aA ^ lnT„ - E(lnT„,) > a A or lnT„ - E(lnT„) < -aX 

^ e^'^^-jC > e'^^ or e^°^"/<^ < e""^- 



Now, fix e > and choose Aq = -^2 In-. Then, by Lemma 0, 

Pr (r„ < Ce~"^ or r„ > Ce"^) < e. 

Theorem |(1) says that if (n + 1)X < Ce-'^^", then E((T„ - {n + l)X)+) is close to 
E(T„ — (n + 1)X). Otherwise, Theorem ^(2) says that Var((T„ — (n + 1)X)+) is bounded. 

Theorem 3 Let C = e'^C'^^"). 

(i) If {n + 1)X < Ce"'"^", then 

|E((T„ - (n + 1)X)+) - E(T„ - (n + 1)X)| < e{n + l)X. 

(ii) If{n + l)X> Ce~''^\ then 

Var((T„ -{n + 1)A)+) <{n + l)^A^e'"^° | ^ 

Ao — 2cr 

Proof: Let 0()f:) denote the probability density function of T^. 
Statement 1. Note first that 

poo poo 

E((T„-(n+l)A)+) = / {t-{n + l)Xy<p{t)dt = / (t - + l)X)0(t)dt 

Jo J{n+l)X 

»(n+l)X 



E(T„ - (n + 1)X) - [ 
Jo 



{t-{n + l)X)(f){t)dt. 
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|E((T„ - (n + 1)X)+) - E(r„ - (n + 1)X) | 



(n+l)X 



< (n + l)X 



(n+l)X 



(p{t)dt 



= (n + l)XPr (T„ < (n+ 1)X) < (n + 1)X Pr (T„ < Ce""^") 
<e{n + l)X, 

where the second-last inequahty follows from the assumption that {n + 1)X < Ce'"'^'^ 
and the last inequality follows from Lemma |^ and our choice of Aq. 



Statement 2. Since (n + 1)X > Ce^'^^o, 
Var((T„ - + 1)X)+) < e(((T„ - (n + 1)X)+)' 

/■CXD 

= 2 / (t - (n + 1)X) Pr(T„ > t)dt 

J(n+l)X 



(t - (n + 



(n+l)X 



< 2 / (t - Ce""^") Pr(T„ > t)dt 



2/ {t-Ce'''^°)FT{Tn>t)dt + 2 (t - Ce~'"^«) Pr(r„ > t)c/t 



Ao 



(substitute t = Ce'^^ in the second term) 



Ao 

g2aAog-A2/2 



Ao - 2a 



(from the inequality e dx < je 2 



2 x^2„4(tAo 



< (n + lyx'e 



2ae 



Ao - 2(T, 

(since (n + l)Xe2'^^" > Ce'^^" and e = 2e"^o/2) 

<(r2+l)2xV'^^«i^^. 
- ^ ^ Ao - 2cT 



5.^. The BoundedMC Algorithm 

We next use these results in our algorithm. One approach would be to estimate C = e^^'^" 
to determine whether we should apply Theorem ^(1) or ^(2). Our algorithm takes a 
more direct approach. We begin by selecting N samples u;^, . . . , tu^ G VL and computing 
T„(u;^), . . . , Tn{ijj^). For 1 < i < N , define the random variable Zi as 




1 li Tn{uj') < [n + l)X , 
otherwise. 



Let Z = Y,i=i Zi. 
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Lemma 4 Let < 6 < 1 be given. With N = 0(log |) trials, the following statements 
hold. 

(i) If§< 2e, then Pr (T„ < {n + 1)X) < 4e with probability 1-S. 

(ii) If§> 2e, then Pr (T„ < {n + 1)X) > e with probability 1 - 5. 

Proof: This follows from the Chernoff bound [|T^, Theorems 4.1 and 4.2]. ■ 

Theorem 5 Let < 5 < 1 be given. With N = 6(logi + ^e^'^^^^J^) trials, the 
following statements hold. 

(i) If§< 2e, then, with probability 1-S, E(r„ - (n + 1)X) estimates E((T„ - (n + 
1)X)'^) with error at most 4e{n + 1)X . Correspondingly, the price of the call is 
estimated with error at most AeX . 

(ii) If > 26, then, with probability 1 — 5, j^J^iLi (^nl"^*) — {n + 1)X)~^ estimates 
E((T„ — {n + 1)X)~^) with standard deviation at most e{n + 1)X. Correspondingly, 
the price of the call is estimated with standard deviation at most eX . 

Proof: The statements are proved as follows. 
Statement 1. With probability 1 — 5, we have 

|E((T„-(n + l)X)+)-E(r„-(n + l)X)| < (n + l)XPr(r„< (n + l)X) 
< Ae{n + 1)X, 

where the first inequality follows from the proof of Theorem ^(1) and the second 
inequality follows from Lemma §(1). Since the price of the call is just ^^E((T„ — (n + 
1)X)^), the error in our estimate for the price is at most AeX. We give a closed-form 
formula for E(T„ — {n + 1)X) = E(T„) — {n + 1)X in Lemma ^ below. 

Statement 2. If ^ > 2e, then by Lemma §(2), with probability 1—5, Pr(T„ < 
(n + 1)X) > e. This implies {n + 1)X > Ce'"^''. Otherwise, if (n + l)X < Ce-"^", we 
would have 

e < Pr(T„ <{n + 1)X) < Pr(r„ < Ce""^"), 

which contradicts our choice of Aq above. Since {n + 1)X > Ce^'^^°, we can 
apply Theorem ^(2). The variance of our estimate ^Ei=i(^n(w*) - {n + 1)X)+ for 
E((T„ - (n + 1)X)+) is ^Var((T„ - (n + 1)X)+). Since 

Var ((r„ -{n + 1)X)+) < (n + ifX^e^'^''^^^^, 

the variance of our estimate is at most + 1)^X^ for X > ■^e^"'^° \^-2a 
standard deviation is e(n+l)X. Since the price of the call is just ^^E((T„ — (n+l)X)~''), 
the expected error in our estimate for the payoff is at most eX. ■ 

We now summarize our algorithm. 
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Algorithm 1 BoundedMC((5, e) 
generate N = Q{\og i + ^e^-'^o^^) paths, w\ . . . , u^; 
let Z be the number of paths u"^ such that T„(ci;*) < {n + 1)X; 
if Z/N < 2e 

return ;^E(r„ - (n + 1)X) = ;^(E(T„) - (n + 1)X); 
else 

return ^Eii(^nM-^)+. 

For the completeness of Theorem ^ we conclude this section with a formula for 
E(T,0 [0. 



Lemma 6 Let 5*0 be tie initial stock price and let r be the interest rate. Then 

{n + l)So, ifr = 0, 
^ — ^ So, ifr>0. 



E(T„) 



Proof: Recall that s{t,i) is the stock price at node T[t,i], u and d are the up- 
tick/downtick factors, and p and q are the risk neutral probabilities. We have 

E(T„)= E s{tJ)-Vi{St = s{tJ)) 



level < t < n node < £ < t 

t=0 £=0 ^ ^ 



5*0 ^^(pu + (1 — (using the binomial theorem) 

i=0 

n 

5*0 ^^(1 + (assuming risk-neutral rate of return) 



{n + l)So, if r = 0, 



(l+r)"+l-l 



So, ifr>0. 



as claimed. 



3.3. Pricing Basket Options using BoundedMC 

The results derived above are applicable for European Asian basket calls as well. 
Suppose we have a basket of m stocks. Let T„ denote the running total of the 
stocks in the basket up to day n, let f2* be the sample space of paths the binomial 
tree corresponding to stock i and let 11* be the corresponding measure on Let 
Q = X ■ ■ ■ X Q"^ and 11 = 11^ x ■ ■ ■ x H™. Lemma § is again applicable, but this time 
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our random variable U = In T„ has mn components. Letting cXmax denote the maximum 
volatility among the m stocks, we can show, as in Lemma |^, that 

Pr {Tn < Ce~''-^^^^^ or T„ > C'e'^--^^^) < 2e"^'/^ 

for any A > 0. The results of the single stock case are now applicable. Given e > 0, 
choose Ao = ^2 Inf, so that Pr (T„ < C'e-'^"--^v^ or T„ > Ce'""-" ^v^) < e. We sample 
paths cu^, . . . in f2 according to 11. Let Z be be number of paths a;* such that 
Tn{'^^) < + 1)-^, where X is the strike price of the basket call. 

Theorem 7 Let < 5 < 1 be given. With N = Gflogi + ^e^"-^-'^ ^ov^ i+2a.nax sV^ ^ 

f V to 5 I Ao-2o-niax V™'^ 

trials, the following statements hold. 

(i) If§< 2e, then, with probability l-S, E(r„ - (n + 1)X) estimates E((r„ - (n + 
1)X)~^) with error at most 4e{n + 1)X . Correspondingly, the price of the call is 
estimated with error at most AeX . 

(a) If > 2e, then, with probability 1 — 6, jj: J2iLi (^n(t^*) — {n + l)X)'^ estimates 
E((T„ — (n + 1)X)"'") with standard deviation at most e{n + 1)X. Correspondingly, 
the price of the call is estimated with standard deviation at most eX . 

Proof: This theorem is analogous to Theorem |^ and can be proven using similar 
techniques. ■ 

4. A Recursive Bucketing-Based Scheme 

The AMO algorithm takes 0{kn'^) time to produce a price estimate in the range 
[P-^,P], where P is the exact price of the call and k is any natural number. As in 
Section ^, this algorithm estimates E((T„ — (n + l)X)+), from which the price of the call, 
K([An — X)^), can be easily estimated. Our recursive scheme is a generalization of a 
variant of the AMO algorithm. Below we first describe this variant, called Bucketed Tree 
Traversal (BTT) and then detail our scheme, called Recursive Bucketed Tree Traversal 
(RecBTT). 

Given binomial tree T of depth n, and numbers t, i, m such that < t < n, 
< i < t, and m < n — t, let Tm^^ be the subtree of depth m rooted at node T[t,i]. 
Given < t < n, let uj\t be the prefix of uj up to level t of T. Note that a;|„ = u. Given 
ipyUJ & Q, we say that ip is an extension of uj\m if, for < t < m, we have ipt = ^t- Given 
another binomial tree U of depth n we say that G VL{U) is isomorphic to G ^{T), 
if, for all Q < t < n, ipt = ^t- 

Like the AMO algorithm, BTT is based on the following simple observation. 
Suppose that the running total Tm{ijj) = Tm{uj\m) of the stock prices on path lu E Q 
exceeds the barrier B = (n+ 1)A. Then, for any extension ip of uj\m, Tn{ip) clIso exceeds 
B and the call will be exercised. If we know the call will be exercised on all extensions 
of u\m, it is easy to compute the payoff of the call on these extensions, as described in 
Lemma ^ 
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As we travel down a path u, once the running total Tm{uj) exceeds B, we can just 
keep track of the running total on extensions ip of uj\m weighted by from which 

the value of the option can be computed. Hence, we only need to individually keep 
track of path prefixes u!\m that have running totals Tm{uj\m) less than B. 

Unfortunately, there may be exponentially many such Lj\m- However, the running 
totals Tm{uj\jn) are in the bounded range [0,B). Rather than trying to keep track of 
each running total individually, we instead group the running totals terminating at each 
node into buckets that subdivide this interval. This introduces some round-off error. 
Suppose we use k buckets to divide [0, B) into equal-length subintervals and we use the 
left endpoint of each interval as the representative value of the running totals contained 
in that bucket. At each step down the tree, when we put a running total into a bucket, 
an error of at most j is introduced. Traveling down n levels of the tree, the total T„(a;) 
of a path u is underestimated by at most ^ and the average An{uj) is underestimated 
by at most 

BTT is detailed in Algorithm 0. At each node v = T[t,i] of T, create k + 1 
buckets to store partial sums of path prefixes terminating at v. There will be k core 
buckets and one overflow bucket. This overflow bucket is the only difference between 
BTT and the AMO algorithm. For < j < k, core bucket bj{v) stores the probability 
mass bj{v). mass of path prefixes that terminate at node v and have running totals in its 
range range(j) = [jj, (j + l)-|). The representative value of partial sums in the bucket is 
denoted by value = j^- The overflow bucket bk{v) stores the probability- weighted 
running total estimates of path prefixes that have estimated running totals exceeding 
B. This quantity is denoted by value. The probability mass of these path prefixes 
is denoted by 6fc(f).mass. BTT iterates through each of the k + 1 buckets of each of 
the nodes in T, for a total runtime of 0{kn'^). 

Algorithm 2 BTT(r, k, B) 

for each node v eT and each bucket bj{v) 
set 6j(w).mass ^ 0; 

take i such that initial price s(T[0, 0]) G range(j); 6j(T[0, 0]).mass <— 1; 7. assume 

s(r[o,o]) < B 

for t = 0, . . . , — 1) °L iterate through each level 

for i = 0, . . . ,t 7o iterate through each node at level t 
let V = 7'[)f:, i]; % shorthand notation for node T[t, z] 
for W e {T[t + l,i],T[t + l,i + 1]} 7. for each child of node v 
let p' G {p, q} be the probability of going from node v to w; 
for bj{v) G {bo{v), . . . ,bk{v)} 7. for each bucket at node v 

let V <— value + s{w); 
let M ^ bj{v). mass x p'; 
ifV <B 

take i such that V G range(£); 

bi{w). mass <— bi{w). mass + M; 
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else y. in overflow bucket 

value ^ fe.(^).massxM^).value+MxV 

6fc(w).mass ^ 6fc(w).mass + M; 
return ^"^g 6jt(n, ■j).mass x (6fc(n, i). value — 5). 

7, return option price estimated from overflow buckets at leaves 

We propose RecBTT, a recursive extension of BTT. Consider some level t in our 
binomial tree T and assume that the weights of all path prefixes terminating at level t 
have been put into the appropriate buckets. BTT uses these weights to compute the 
bucket weights of nodes at level t + 1. In contrast, RecBTT recursively solves the 
problem for subtrees 7^*'*^, < i < t, of some depth m < n — t rooted at node T[t, i].^ 
As each recursive call is complete, RecBTT Merges the bucket weights at the leaves 
of 7^*'*' into the corresponding nodes at level t + m oiT. The advantages of the recursive 
calls are twofold. 

(i) They use finer bucket granularity, resulting in improved accuracy. 

(ii) The results of a single recursive call on a particular subtree T^*'*' are used to 
Estimate the results of other recursive calls to other subtrees Tm''\ where j > i, 
as long as the node prices in Tm''^ are "sufficiently close" to the corresponding node 
prices in 7^*'*'. This improves the runtime, since we do not need to make alH + 1 
of the recursive calls, so there are portions of T that we do not directly traverse. 



4.1. The Merge Procedure 

Consider a recursive call on the subtree Ti = 7^'*"'*°' of depth ni rooted at node 
Vq = T[to,'^o]- A leaf vi = Ti[ni,ii] of 7^ (0 < ?i < ni) corresponds to the node 
V2 = T\to + ni,zo + ^i] of 7". Merge incorporates the bucket weights at vi into the 
bucket weights at V2- Recall that the recursive call on 7j is made with finer bucket 
granularity. Assume we use ki = hiko core buckets instead of just = k. We first 
combine each group of hi buckets at vi into a single bucket, so that we are left with 
ko buckets to merge into V2. When we refer to a bucket below, < j i < ko, we 

mean one of these ko combined buckets. 

The Merge procedure is described in Algorithm ^. Consider first the core buckets. 
Let < jo,ji,j2 < ko denote core bucket indices. Bucket bj^Xvo) contains the mass of 
path prefixes in T terminating at vq whose running total estimates fall into the interval 
range(jo) = [jo^; (jo + 1)^)- Bucket &ji(fi) contains the mass of full paths in 7^ 
terminating at node Vi whose running total estimates fall into the interval range(ji). 
This is equal to the mass of partial paths in T starting at node Vq and terminating 
at node V2 whose running total estimates fall into the interval range(ji). Merging 7i 
into T involves merging each leaf Vi of 7^ into the corresponding node f 2 of T. Once 

+ Actually, the trees on which we recursively solve the problem are not exactly the 7^*'*'. Each is 
identical to the respective Tm^\ except the price at the root is changed from s(T[t, i]) to 0. The reason 
for this is explained in Remark |^. 
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the merging procedure is done, fc^-j (1)2 ). mass is updated to contain the weight of path 
prefixes in T passing through node Vq and terminating at node V2 that have running 
total estimates in the interval range(j2)- The overflow buckets are handled similarly. 

Algorithm 3 MERGE(r, Ti = Xt'""^) 
let vq T[to,io]; 

for ii = 0, . . . , rii 7. for each leaf of 7i 

let I'l ^ Ti[ni,ii], V2 T[to + ni,io + ii] ; 

for jo — 0, . . . , ko 7. buckets in vq 
for ji = 0, . . . , ko 7. buckets in vi 

let V <— 6jo(wo)-value + value; 
let M <— 6j(,(t'o).mass x 6j^(t'i).mass; 
iiV <B 

take j2 such that V e range (J2); 
(1)2 ). mass <— ('^2) -mass + M; 

else 



6/5(^2). value 



^fe(''2)-iiiassxf)j.(i;2).value+Mxy _ 



fefc(r)2).mass+M 

6fe(v2).mass <— 6jk(v2).mass + M . 



Remark 1 Notice that the price at Vq is counted twice: once in the path prefix from the 
root T[0, 0] to Vq and once in the partial path between Vq and V2. To address this issue, 
when we recursively solve the problem on the subtree Ti, we set the price s(Ti[0,0]) at 
the root to be 0, ensuring that this price is counted once. This modification does not 
change our algorithms. 



Lemma 8 For an arbitrary node v, let E{v) he the maximum amount by which running 
totals terminating at node v are underestimated by the bucket values. Using the Merge 
algorithm, we have E{v2) < E{vq) + E{vi) + j. 

Proof: Consider first the core buckets. Bucket bjg{vQ) covers running totals in the 
range {joj, (jo + l)f ) and 6^1(^1) covers running totals in the range [(ji)f , (ji + l)f ). 
The concatenation of paths in these buckets covers running totals in the range 
[(jo+Ji)f , (io+ii + 2)f ). Assuming jo+ji < k, we put these paths in bucket bj^+j^{v2), 
and the running totals are underestimated by the original error of E{vo) + E{vi), plus 
an additional ^, since some partial paths that fall in the range of bucket 6jg+j^+i(f2) are 
actually put in 6^0+^1(^2). If Jo+Ji ^ k, then these paths are put in the overflow bucket. 
In this case, no additional error is introduced. The original error of E{vq) + E{vi) is 
carried over. Similar arguments can be made when we are merging overflow buckets. ■ 



Lemma 9 Merge can be made to run in 0{niklogk) time. 
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Proof: In the implementation of Algorithm ^ MERGE runs in 0{nik'^) time. 
However, the core buckets can be merged with a faster technique. Let ao{x) = 
E,o~=o^io(^^o).mass-x^« andai(x) = E,-, =0 ^ii ('Wi)-niass-X''^ be polynomial representations 
of the bucket masses at nodes Vq and Vi. Let 02 (x) = ao{x) ■ ai{x). This product can be 
computed in 0{klogk) time with the Fast Fourier Transform (FFT) |^. For < j 2 < k, 
the coefficient of the x^'^ term of 02 (x) is the probability mass that should be added to 
bucket hj^{y2). For j2 > k, this coefficient should be added to the overflow bucket hkiy-z). 
So the core buckets can be merged in 0{k\ogk) time per node. Merging the overflow 
buckets takes 0{k) time, so the total runtime per node is 0{k\ogk). ■ 



4-2. The Estimate Procedure 

Let Ti = Tm°'^^^ and T2 = Tm^'^'^^ be two subtrees of T, where ^2 > ii We now describe 
the Estimate(7^, Ti) procedure, which estimates the weights in the leaf buckets of 7^ 
from the weights in the leaf buckets of Ti. This saves us the work of recursively solving 
the problem on 7^. 

Estimate is described in Algorithm ^. It uses the following fact. Given any node 
Vi = Ti[t,i] in Ti, let V2 = T2[t,i] be the corresponding node in 7^. Notice that there 
is a constant a > 1 such that for all (fi,f2) pairs, 5(^2) = as(fi). Hence, for any path 
ip G ^2(72), we have Tm{ip) = aTm{uj), where u G fl(7i) is isomorphic to ip. 

Algorithm 4 Estimate(T2 = Ti*"'*'', Ti = T^''''^) 

for i = 0, . . . ,m 7, go through the leaf buckets of Ti and T2 

let Vi ^ Ti[m,i], V2 ^ 7^[m,i]; 

for ji = 0, . . . , 7o go through each bucket at vi 

let V <— afojj value; 
let M ^ 6j^(t>i).mass; 
if < 5 

take j2 such that V G range(j2); 

bj2{v2)-iaa.ss ^ 6^2(^2). mass + M; 
else 

h (1, ^ iralne ^ fefc (f 2 ) .mass X &fc (i>2 ) . value+M X V . 
Ofcl^^2j-VaiUe < 6fc(t,2).mass+M ' 

6fc(f2).mass <— bk{v2) -Taass + M; 



Lemma 10 Suppose that a < 2 and assume that the total path sums in 7i are 
underestimated by our bucketing scheme by at most E. Estimate underestimates 
the total path sums in 7^ by at most 2E + 2 

Proof: Take any path u G ^l{Ti) and let Um{'^) be the estimate for Tra{oj) made by the 
bucketing scheme. By assumption, Tmiyj) — Umiyj) < E. Let ip G i^(7^) be isomorphic 
to oj. Using aUmiyj) as an estimate for Tmiip) = aTm{uj), we underestimate T^i^ip) by 
at most aE < 2E. This accounts for the flrst error term. 
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Now, consider core bucket ^^^(fi) of some leaf Vi of The range of running totals 
covered by this bucket is range (ji) = [jif , (ji + l)-f ). When we map this bucket to 72, it 
covers path sums in the range [«ji-f , + l)"f )> which falls into at most three buckets 
in the corresponding node V2 of 7^, since a < 2. Call these buckets hj^{y2), 6^2+1 ("^2), 
and hj^+2{v2) (assume these are all core buckets). Putting the entire mass into hj^{v2) 
causes the portion of partial sums that fall into bucket 6^2+2(^2) to be underestimated 
by additional 2^, which accounts for the second error term. 

Similar errors occur if the range [ajif + 1)^) spills into the overflow bucket 
bkiv2) and when mapping the overflow bucket hkivi) in t>i to 7^. ■ 

Lemma 11 Suppose that we would like to determine the leaf bucket weights of the 
subtrees Tm^\ where < i < t. We need only call RecBTT at most once for every 
0(^) subtrees, and use the Estimate procedure to estimate the leaf bucket weights 
of the other subtrees with bounded error. 

Proof: We begin by recursively solving the problem on the subtree Tm^K Notice that 
s{T[t,i]) = M^*s(T[t, 0]). By Lemma as long as m^* < 2, we can use Estimate to 
compute the leaf bucket weights of 7^*'*^ , without introducing significant error. Recalling 
that we take u = e"'^^, where a is the volatility of the stock, we see that we can use 



Estimate (7^*'*', 7^*'°') for i as large as ^^'^^ = 6(^) while keeping the error bounded. 



4-3. Error and Runtime Analysis 

We now derive recursive expressions for the error and runtime of RecBTT. Suppose 
there are Hq trading periods and we use a total of ko buckets per node. For i > 0, let 
Ui be the number of trading periods and ki the number of buckets per node we use in 
the i-th subproblem (at the i-th level into the recursion). Let Ei = E{ni,ki) be the 
error bound made in the i-th subproblem. The error is the amount by which RecBTT 
underestimates the running totals. Let Tj = T(nj,/cj) be the runtime of RecBTT on 
the i-th subproblem, including all recursive calls. We recursively call RecBTT on 
binomial trees of decreasing depth (decreasing ni) with an increasing number of buckets 
(increasing ki) until we reach a tree of some small depth. 

Lemma 12 When i recursive calls are made, the error made by RecBTT is 



where Ei+i is the error in solving the {i + l)-st subproblem. 

Proof: Consider the error in the i-th subproblem, Ei = E{ni,ki). Each recursive 
call has an error of -Ej+i. When the fcj+i smaller granularity buckets in the {i + l)-st 
problem are grouped into the ki < /cj+i buckets of the i-th problem, an error of one 



This proves the lemma. 




j=0 
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bucket size, ^, is introduced. A call to Estimate underestimates the price sums by 
at most 2(Ei^i + + 2^ = 4^ + 2£',;_|_i, by Lemma 0. This error is made at most 
times, once at each of levels n^+i, 2nj+i, . . . , ;^^nj+i. At each of these levels, an 
additional error of ^ is introduced by the Merge procedure, as described in Lemma |. 
The total error in the i-th subproblem is 

Ei = 1 -c/i+i- 

Unraveling i levels of the recursion on the original problem yields the result. ■ 

For the runtime of RecBTT, we must analyze how many times the {i + l)-st 
subproblem needs to be solved when solving the i-th subproblem. At each of levels 
0,nj+i, . . . , — in the binomial subtree of the i-th subproblem, suppose the 

{i + l)-st subproblem is called at most q times. 

Lemma 13 When i recursive calls are made, the runtime of RecBTT is 
To = o(j2 non.k, log(fc,) j] c)\ + ^^^^D^T.+i, 

Vj=0 k=0 J 

where Tj+i is the time it takes to solve the {i + l)-st subproblem. 



Proof: Consider the runtime of the i-th subproblem, Tj = T{ni,ki). Each 
recursive call takes time Tj+i and at most Cj recursive calls are made at each of levels 
0,nj+i, . . . , — The leaf bucket weights of the other at most Ui subtrees 

rooted at the nodes at each of these levels can be determined using Estimate, which 
takes C(fcinj+i) time. For the Merge procedure, at each of levels t, each of the at 

most Ui nodes at level t must be merged with nj+i leaves from the subtrees Tn^fl to get 
the bucket weights at level t + nj+i. The merging of each node takes 0{ki log /cj) time. 
The total runtime of the i-th subproblem is 

Tl' TIC' 

Ti = — —{ciTi+i + 0{nikini+i) + 0{nini+iki\ogki)) = -^Ti+i + 0{n^kilogki). 
Unraveling i levels of the recursion on the original problem yields the lemma. ■ 

Theorem 14 Given integer R > 2, let 'j = j^, and for i > 0, let ui = (^)^/^"*'^ 
and ki = A'^koC^Y'^ , wiiere a is the volatility of the stock. RecBTT underestimates 
E((T„-(n + l)X)+) byatmostC(:^^^2^!!!£ll!Z) and takes time O {2^ /^nlko{^ + hg!^)) . 

Proof: After K = ^ levels of recursion, we have = (^)° = 1- The problem 
consists of a single node and can be solved in constant time Tk with an error of 
Elf = -r- = — /o . ; From Lemma O, the error is 



hBun > 1 Ek = hBun > . ^ , ,^ ,^ h 



1/2 



o\- ^ 

fco 
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For the runtime analysis, by Lemma |Tl|, note that Cq = = a^Jn. For i > 0, since 



rij < Cj = 1. Combining this fact with Lemma |T3| yields the runtime 

noUiCokilogki 1 + 0(?t,oCo) 



, i=l 



nlko logko + OiY. ^o4'^o(^ + log ^ 
\i=i 

0(2V^n2fco(i + log^)). 



Corollary 15 Given integer R> 2, let j = and choose Ui and ki as in Theorem 
In time 0{2^/^nlko{- + log^)), RecBTT returns an option price in the range 

[P ^^~~k^^ ' -^]' ^here P is the exact price of the option. 

Proof: This follows from Theorem O and the fact that 



E((A„ - (n + 1)X)+) = ^E((T„ -{n + l)Xy 

n + 1 



Remark 2 By comparison, the AMO algorithm |l|] has error (9(noX/A;o) and runtime 
O{konl). 

4.4- Customization of the RecBTT Scheme 

A useful feature of the RecBTT scheme is that the number of recursive calls, the 
size of each recursive call, and the pricing method by which we solve the base case 
of the recursion (the final recursive call) can be custom tailored for the application. 
The base case of the recursion can be solved with any option pricing scheme, including 
BoundedMC, other variants of the Monte Carlo algorithm, the AMO algorithm, or 
exhaustive path traversal. For example, consider a 3-month (90-day) Asian option with 
the average computed on a daily basis (n = 90). Pricing this option exactly requires 
traversing 2^° paths, which is infeasible. Using our RecBTT scheme, since y/n < 10, a 
single recursive call can reduce this problem to several subproblems each with at most 
2^° paths. 

In practice, larger values of n appear in several applications. Long-term options 



contracts, LEAPS |jT3[, are negotiated for exercise dates several years in advance. 
Companies also offer stock options to employees over periods of up to five years or 
more. Finally, in Istanbul options contracts rather than computing the average 
based on daily prices, the average is computed over prices taken at higher frequencies. 
In each of these cases, n is sufficiently large that several recursive calls are required to 
reduce the problem to a manageable size. 
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5. Pricing Basket Options 

The bucketed tree structure created by BTT and RecBTT can be used to price various 
kinds of European basket options as well. Here, we describe here how to price European 
Asian basket options. A basket option |T3| is composed of m stocks, zi, . . . ,Zm- For each 
stock Zi, we construct a binomial tree (according to the respective stock's volatility), as 
described in Section |^. For 1 < t < n and 1 < i < m, let SI be the random variable 
denoting the price of stock Zi on day t and define St = J21Li SI to be the random 
variable denoting the total of the stock prices on day t. Recall that the payoff of a 
European basket call with strike price X is ]E((S'„ — X)~^). Letting An = Yl^=o'St 
be the average total stock price, the payoff of a European Asian basket call with strike 
price X is ]E((y4„ — X)'^). There is additional complexity with pricing basket options 
that does not appear in their single-stock counterparts: the number of paths that the 
total basket price can follow is exponential, not only in the number of trading periods n 
but also in that of stocks m. Basket options are usually priced using traditional Monte 
Carlo methods. The scheme we describe here is the first polynomial time, in both the 
number of stocks and trading periods, pricing scheme for any kind of basket option with 
provably small error bounds. 

Our European Asian basket call pricing algorithm, BasketBTT, is described in 
Algoirthm^. Let B = {n+l)X, where X is the strike price of the basket option. For each 
stock Zi, 1 < i < m, use RecBTT to construct the bucketed binomial tree structure T* 
described in Section ^, this time using B as the barrier; should the running total of any Zi 
exceed B, the basket option will always be exercised, regardless of what the other stocks 
do. For each stock Zi, we construct k + 1 superbuckets l3j, < j < k, where /3j is the 
combination of buckets bj{v) for all leaves f G T*. For the core buckets Pj, < j < k, let 
value = jj and /5j.mass = ^"^q 6j(T*[n, £]).mass, where this summation ranges over 
all leaves T*[n, £] of T\ For the overflow bucket PI, let /3^.mass = X]"=o £]).mass 
and 

ELo^fc(^1^'^]) -value X bj{T'[n,i]). mass 



Pl-valne 



/SL.mass 



Handling overflow superbuckets If the running total of a stock Zi reaches the overflow 
superbucket PI, the option will be exercised regardless of what the other stocks do. 
Given this, we can determine the value of the option exactly, since 

E((T„ -{n + 1)X)+) = E(T„ - (n + 1)X) = E(T„) - {n + 1)X 

= value + ^ E(T^') - + 1)^, 

where is the random variable denoting the running total of stock Zi/ up to day n. 
E(T^') can be computed exactly using Lemma ^ 

Handling core superbuckets Consider now the core superbuckets Pj, < j < k. Let 
fi{x) = ^^~Q/9*.mass ■ be the polynomial representation of the core bucket masses 
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of stock Zi and let /(x) = YliLifii^)- '^^^^ product can be computed efficiently, as 
described in Lemma |16|. Notice that f{x) has the form 

^From the definition of f{x), observe that bj is just the probability that the sum (over 
all stocks Zi) of running totals from the core buckets falls in the range [j-^, (j + l)f )• 
That is, 

XjT:G[jf,(j + l)f) |T:<5for alH 
Hence, the contribution to the option price from the core buckets can be estimated by 

m{k—l) 

J2 b,{jf-{n + l)X). 

j=k 

Pricing the option Combining the above results for the overfiow and the core 
superbuckets, we see that IE((T„ — {n + 1)X)~^) can be estimated by 

m m(fc-l) 

^/3^mass(/3^value + E.v.E(^n)-(^ + lW+ E b,{jf-{n + l)X). 

1=1 j=k 



Algorithm 5 BasketBTT(zi, ...,Zm,B = {n + 1)X) 

for i = 1, . . . , n 7, for each stock Zi 
compute lE(T^) using Lemma ^ 

run RecBTT on stock Zi with barrier B and k buckets; 
for J = 0, . . . , — 1 y, construct core superbuckets 



B. 

k ' 



value j 
/Sj.mass <— ^"^Q 6j(T*[r;,, £]).mass; 

mass ^"^Q 6j(T*[n, £]). mass; 7. construct overflow superbucket 

value ^ ^. .rlass XlLo ^fc(^'[^, ^]). value x 6j(T*[n, £]). mass; 
let 4- ^^^Q/3*.mass x^; 

rm 



compute /(x) ^ Ili^i /«(^) described in Lemma |T^; 

let f{x) ^ boX° + biX^ H h 6m(fc-l)a;™^''"^''; % for some 6o,---,&m(fc-i) 

return ^ /5^.mass(/3^value + E.v. ^(^^V (^ + 1)^) + Erir'^ ^.O'f - 1^^ + 1)^); 



Lemma 16 Let (Tmin be the minimum volatility among the stocks in the basket. The 
runtime of BasketBTT is 

(9(mT((Tmin ,n,k) + mk log m log k + mk log^ m), 

where T((Tmin , n, k) is the runtime of RecBTT on a binomial tree of size n when k 
buckets are used and the volatility of the underlying stock is cTmin ■ 
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Proof: Computing the leaf bucket weights for each stock Zi (using RecBTT) takes 
0{m ■ T(crinin ; k)) time. The runtime of the rest of the computation is dominated by 
the time to compute YYiLi fi{^)- We now describe an efficient way to do this. Assume 
that m = 2*", for some r. The general case is handled similarly. We conduct the 
multiplication by successively multiplying consecutive polynomials together logm = r 
times until we are left with a single polynomial. For 1 < z < m, let fi{x) = fi{x). 
For 1 < j < r and 1 < « < let fi{x) = f2i^i{x)f2i~'^{x)- The answer that we 
are looking for is just YULi fi{^) = fii^)- At stage j, we multiply together ^ pairs 
of polynomials, each of degree 2^~^k. Using FFT, this takes 0{^2^~^k\og{2^~^k)) = 
0{mklog{2^~^k)) time. The total runtime to compute the product, over all stages, is 
0{mk{r\ogk + X]j=i log2-'"^)) = (9(mA;(logmlogA; + log^m)), from which the claimed 
result follows. ■ 

Our definition of the error made by BasketBTT is symmetric to the definition of 
error made by RecBTT; i.e., it is the maximum amount by which BasketBTT can 
underestimate J2t=o ^ti^^)^ paths cj* G T\ 

Lemma 17 Let (Xmax be the maximum volatility among the stocks in the basket. The 
error made by BasketBTT is at most mi?((Jmax ,n,k), where E{a^a_x , n, k) is the error 
made by RecBTT on a single stock with volatility amax ■ 

Proof: For 1 < i < m, let Ui be any path down the binomial tree corresponding 
to stock Zi and let T^{u)i) be the total price down Ui. When we run RecBTT on 
Zi, each T^{uJi) is underestimated by at most -E(crmax ,n,k). Hence for any uji, . . . , ujmi 
Y^i'^ni'^i) is underestimated by at most mE{auia.x,n,k), as claimed. ■ 

Theorem 18 Given n, m, k, R > 2, 'j = and o"max, if we apply RecBTT 

as described in Theorem [IJ to construct the bucketed binomial tree for each stock, 
BasketBTT has an error of (^(m ^"^'^^'^'^'"'^^^ ) and runs in time 

0{2^^"'n'^mk(- + log -^^) + mA; logm log k + mk log^ m). 



Proof: This follows directly from Lemmas 16 and O and Theorem M. 



Corollary 19 Given n, m, k, R > 2, = j^, a^-m and (Jmax, BasketBTT 
underestimates the price of a European Asian basket call by at most (9(m^^^^-^^^p^s^) 
and runs in time 

0{2^/'^n^mk{- + log -^^) + m/c logm log k + mk log^ m). 



mm 



'7 

Proof: This follows from Theorem [T^, which gives the amount by which the total 
prices of the stocks may be underestimated. The amount by which the sum of their 
daily averages is underestimated is smaller by a factor of ■ 
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6. Further Research 



This paper focuses on theoretical aspects of BoundedMC and RecBTT. In order 
to optimize the performance of RecBTT in practice, we must consider three issues 
in its implementation: (1) The number of recursive calls and the size of the recursive 
subproblems will depend on the number of periods (see Section [4.4|) . (2) For options 
with a high number of periods, nontrivial data structures may be required. (3) We have 
a choice of several algorithms to use at the base-case of our recursion (see Section [4.4|) . 
The implementation of our algorithms and empirical comparisons between them and 
others will be reported in a subsequent paper. 

This paper opens up several directions for further theoretical research. An 
immediate open problem is whether similar Monte Carlo and bucketing techniques can 
be used for pricing American Asian single-stock and basket options. 

With regards to our algorithms, it would be an important result if the runtime of 
BoundedMC on basket options can be bounded by a polynomial in the number of 
stocks m; or if the recursive calls in RecBTT can be structured differently to yield 
improved runtime or error bounds. 

Finally, it would be interesting to see whether our Monte Carlo analysis or recursive 
bucketing techniques can be used to price other types of options or applied to similar 
problems. 
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